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Abstract 

The exclusive photoproduction reaction 7 p — > T p has been studied with the 
ZEUS experiment in ep collisions at HERA using an integrated luminosity of 
468 pb" 1 . The measurement covers the kinematic range 60 < W < 220 GeV and 
Q 2 < 1 GeV 2 , where W is the photon-proton centre-of-mass energy and Q 2 is the 
photon virtuality. These results, which represent the analysis of the full ZEUS 
data sample for this channel, are compared to predictions based on perturbative 
QCD. 
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1 Introduction 



Exclusive photoproduction of heavy vector mesons, J/ip and T, has previously been stud- 
ied at HERA [1-4]. The process 7 p — > V p, with V = J/ip,T can be described by 
perturbative QCD (pQCD), since the relatively high masses of the charm and bottom 
quarks provide sufficiently hard scales. The process is assumed to exhibit threefold fac- 
torisation [5-8]: the photon fluctuates into a qq pair; the pair interacts with the proton; 
and finally the heavy meson is formed in the final state. At leading order (LO), the 
interaction of the qq pair with the proton proceeds via the exchange of two gluons in 
a colour- singlet state. Thus the cross section is proportional to the square of the gluon 
density in the proton. The rise of the gluon density with decreasing fractional momentum, 
x, leads to the prediction of a cross section rapidly rising as a function of the photon- 
proton centre-of-mass energy, W, where the relevant x region accessible in heavy-quark 
production at HERA is 10~ 4 < x < 10~ 2 . The rise of the cross section, a, with W can be 
expressed as a oc W 5 , where a value of 5 ~ 0.7—0.8 has been measured for the J/ip [1,4]. 
A value of 1.7 is predicted for T(1S) production in LO [5]. 

Prior to this analysis, ZEUS and HI measured the T photoproduction cross section for one 
value of W. The increased statistics of the data used in this study allows the investigation 
of the dependence of the production cross section on the energy W. The data cover the 
kinematic range 60 < W < 220 GeV and the results are obtained using the decay 
channel. In this measurement the three upsilon states T(1S), T(2S) and T (3S) (denoted 
Tj, 2=1,2,3, respectively) are not resolved. Hence the sum of the cross sections multiplied 
by the corresponding decay branching ratios to muons, ^\cr 7P_>TiP ■ £>,j, was measured 
as a function of W for two intervals: 60 < W < 130 GeV and 130 < W < 220 GeV. 
The sample under study represents more than a ten-fold increase in integrated luminosity 
compared to the previous ZEUS publication [2]. 

2 Experimental set-up 

In 1998-2007 (1996-1997), HERA provided electron^ beams of energy E e = 27.5 GeV 
and proton beams of energy E p = 920 (820) GeV, resulting in centre-of-mass energies of 
y/s = 318 (300) GeV, corresponding to integrated luminosities of 430 (38) pb _1 . 

A detailed description of the ZEUS detector can be found elsewhere [9]. A brief outline 
of the components that are most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central 
1 Electrons and positrons are both referred to as electrons in this paper. 
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tracking detector (CTD) [10] and, for the data taken after 2001, also in the microvertex 
detector (MVD) [11]. These components operated in a magnetic field of 1.43 T provided 
by a thin superconducting solenoid. The CTD consisted of 72 cylindrical drift chamber 
layers, organised in nine superlayers covering the polar-angk^l region 15° < 9 < 164°. The 
MVD provided polar angle coverage from 7° to 150°. 

The high-resolution uranium-scintillator calorimeter (CAL) [12] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/ \f~E for electrons 
and a(E)/E = 0.35/ VE for hadrons (E in GeV). 

The muon system consisted of barrel, rear (B/RMUON) [13] and forward (FMUON) [9] 
tracking detectors. The B/RMUON consisted of limited-streamer (LS) tube chambers 
placed behind the BCAL (RCAL), both inside and outside the magnetised iron yoke 
surrounding the CAL. The barrel and rear muon chambers covered polar angles from 34° 
to 135° and from 135° to 171°, respectively. The FMUON consisted of six planes of LS 
tubes and four planes of drift chambers covering the angular region from 5° to 32°. The 
muon system exploited the magnetic field of the iron yoke and, in the forward direction, 
of two iron toroids magnetised to 1.6 T to provide an independent measurement of the 
muon momentum. 

The iron yoke surrounding the CAL was instrumented with proportional drift chambers 
to form the Backing Calorimeter (BAG) [14,15]. The BAC provided position informa- 
tion for muon reconstruction with a resolution of few mm for two of its coordinates and 
approximately 15 cm for the third, covering the full angular range. 

The luminosity was measured using the Bethe- Heftier reaction ep — > p with the lu- 
minosity detector which consisted of independent lead-scintillator calorimeter [16] and 
magnetic spectrometer [17] systems. 



2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The polar angle, 9, is 
measured with respect to the proton beam direction. 
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3 Kinematics 



The four-momenta of the incoming electron and proton, and the scattered electron and 
proton are denoted by k,p, k! and p', respectively. The exclusive reaction 

ep — > eTp — > e[i + [i~p (1) 

at a given s = (k + p) 2 for electrons and protons, is described by the following variables: 

• Q 2 = —q 2 = —(k—k f ) 2 , the negative four-momentum squared of the exchanged photon; 

• V = (<? ' p)/{k ' P)i the fraction of the electron energy transferred to the hadronic final 
state in the rest frame of the initial-state proton; 

• W 2 = (q + p) 2 = —Q 2 + 2y(k ■ p) + M 2 , the centre-of-mass energy squared of the 
photon-proton system, where M p is the proton mass; 

• M M + M -, the invariant mass of the pair. 

Selected events (see Section 4) are restricted to Q 2 values from the kinematic minimum, 
Qmm = M 2 y 2 /(1 — y) « 1CT 9 GeV 2 , for y = 0.2, where M e is the electron mass, to a 
value at which the scattered electron starts to be observed in the CAL, Q^ax ~ 1 GeV 2 . 
The median Q 2 value is 10~ 3 GeV 2 and the photon-proton centre-of-mass energy can be 
expressed as 

W 2 ^AE p E e y^2E p (E-p z ) 1 (2) 

where (E — pz) is the difference between the energy and the longitudinal momentum of 
the pair. 

4 Event selection 

Exclusive events in photoproduction were selected using dedicated triggers and 

offline selection cuts. At the trigger level, at least one CTD track matched with a 
F/B/RMUON deposit was required. The offline selection was [18]: 

• CAL timing and vertex position consistent with a nominal ep interaction; 

• two oppositely charged tracks matched to the vertex and no other track in the central 
tracking system; 

• at least one track identified as a muon according to a procedure which uses information 
from B/RMUON, FMUON or BAC, whenever available in a given event [19]; if not 
explicitly identified as a muon, the second track had to be consistent with a minimum- 
ionising particle; 
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• tracks with hits in at least 5 CTD superlayers, to ensure a good momentum resolution; 

• \Vi ~ V2 1 < 1-5, where r/i is the pseudorapidit}{§ of a given track, to reduce the influence 
of the purely electromagnetic Bethe- Heftier background; 

• transverse momentum of a track pt > 1.5 GeV; 

• a cut \tt — 0i — 2 \ > 0.1, where 9i is the polar angle of a given track, to further reject 
cosmic-ray events; 

• invariant mass M M + M - in the range between 5 and 15 GeV; 

• the energy of each CAL cluster not associated to any of the final-state muons to be 
less than 0.5 GeV. This threshold was set to be above the noise level of the CAL. It 
implicitly selected exclusive events with an effective cut Q 2 < 1 GeV 2 ; 

• to suppress the contamination from proton-dissociative events, ep — > eTY, the sum of 
the energy in the FCAL surrounding the beam hole had to be smaller than 1 GeV [18]. 
This corresponds to an effective cut on the mass My of the dissociated system, 
M y < 4 GeV. 

The events were selected in the kinematic range 60 < W < 220 GeV. 



5 Monte Carlo simulation 

The detector and trigger acceptance and the effects due to imperfect detector response 
were determined using samples of Monte Carlo (MC) events. Vector-meson production 
was simulated using the DIFFVM 2.0 generator [20] which is based on Regge phenomenol- 
ogy and the Vector Dominance Model in which the photon fluctuates into a virtual vector 
meson which interacts diffractively with the proton via Pomeron exchange. For exclu- 
sive vector-meson production, s-channel helicity conservation (SCHC) was assumed. An 
exponential dependence ~ e bt was assumed for the differential cross section in the four- 
momentum-transfer squared distribution at the proton vertex, t, with a slope parame- 
ter b = 4.5 GeV -2 , consistent with the value obtained for exclusive J/ip electroproduc- 
tion [1,4]. The W dependence of the ■yp—^Tp cross section was parameterised as cx IV s , 
with_| 5 = 1.2. Radiative corrections are of the order of 1 % [21] and are not included 
in the simulation. The DIFFVM generator supplemented by the JETSET 7.3 package 
[22] was used to estimate the value of My below which background events from proton 
dissociation, yp — > TV, do not leave deposits in FCAL. 

3 Pseudorapidity is defined as rj = — In (tan |). 

4 This value was obtained from an iterative process of the analysis. 
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The background, consisting of Bethe-Heitler dimuon events, a purely electromagnetic 
process, was simulated using the GRAPE vl.lk [23] MC program. All MC events were put 
through the simulation of the ZEUS detector based on the GEANT [24] program versions 
3.13 (1996-2000) and 3.21 (2003-2007) and were analysed with the same reconstruction 
and offline procedures as the data. In addition, corrections [18, 19] of the muon detector 
efficiencies determined from a data set consisting of isolated J/ip and Bethe-Heitler events 
were applied. 

6 Signal extraction 

The distributions of the invariant mass for the selected events are presented in 

Fig. Q] for the two W intervals, 60 < W < 130 GeV and 130 < W < 220 GeV, and for the 
entire sample. A clear signal is seen around 9.5 GeV. 

The mass resolution in the T resonance region is approximately 0.2 GeV and does not al- 
low the Tj, z=l,2,3, states to be resolved. The number of signal events in each W interval 
was determined using the following method. The MC of the Bethe-Heitler process was 
normalised to the data in a mass window not containing resonances (5 to 9 GeV, 10.7 to 
15 GeV). The spectrum outside the resonance regions was reproduced and the distribu- 
tion was used to obtain the signal by subtracting the background under the resonances. 
It was assumed that the above three states are produced at HERA in the proportions 
0.73:0.19:0.08, as measured in pp collisions at the Tevatron [25]. The extracted number 
of signal events in the signal region 9-10.7 GeV are tabulated in Tabled] for the different 
W ranges. The DIFFVM generator was used to obtain the combined shape of the three 
Tj states (see Fig. 1). The DIFFVM MC was normalised to the observed number of 
signal events. Using the above-mentioned ratios, the number of T(1S) candidates was 
calculated and is also presented in Table [TJ 

The procedure to determine the fraction of proton-dissociative events in the final sam- 
ple, /pdiss, has been described in detail elsewhere [1]. Due to the small rate of proton- 
dissociative T production, it was not possible to use the analysed data themselves to 
evaluate this fraction even on the full data sample. However, this contribution is ex- 
pected to be similar in all diffractive vector meson production processes [26]. Therefore, 
diffractive J/ip meson production [1] was used to evaluate this correction. The value 
assumed in this paper was / p di SS = 0.25. 

The effective photon flux calculation [27] takes into account the transverse and longi- 
tudinal flux factors [20]. Thus, the effective flux calculation depends on the assumed 
parameterisation W s of cr{^p). The uncertainty on the W dependence is taken into ac- 
count clS db systematic uncertainty. 
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7 Systematic uncertainties 



The following sources of systematic uncertainty were considered [18] (numbers for the full 
W range are given): 

• signal extraction method: the background was normalised separately to the left and 
to the right sides of the mass spectrum outside of the signal region: 1 9 %; 

• the uncertainty of the CTD tracking and muon-chamber performance in the trigger, 
and the subsequent muon reconstruction in the offline analysis: ±9%; 

• /pdiss wa s varied between 0.2 and 0.3 to account for uncertainties in the experimental 
conditions and for the fact that it was evaluated from J/ifj production [1], resulting in 
an uncertainty of ±7%; 

• the uncertainty related to a variation of the treatment of the calorimeter noise: —4%; 

• the variation of 5 = 1.2 ± 0.5 results in an uncertainty of +3.2, —2.2% on the ep cross 
section and —3.9%, +1.6% on the 7p cross section; 

• variation of the slope parameter b = 4.5 ± 0.5 GeV~ 2 gives a negligible effect on the 
acceptance. 

In addition, the following uncertainties related to the T(1S) extraction were also consid- 
ered: 

• the T(1S) fraction in the signal was varied within the total uncertainties quoted by 
the CDF collaboration [25]: +5%, -6%; 

• the T(1S) decay branching ratio: ±2.4% [28]. 

The total systematic uncertainty was determined by adding the individual contributions 
in quadrature. The values for the different W ranges are given in Table [TJ 

The uncertainty of the luminosity determination is 2.6 % and is not included in the result. 

8 Results 

All ep cross sections are given at yfs = 318 GeV. 

For each W bin, the sum of ep cross sections multiplied by the corresponding decay 
branching ratios to muons was evaluated according to the formula 




(3) 
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where N sig is the number of signal events in the signal mass region, A is the overall accep- 
tance, Bi is the decay branching ratio into and C is the corresponding integrated 
luminosity. 

The ep cross section for T(1S) production was calculated according to 

B\ 

i 

where f r{is) = 0.73±g;g| [25] and Bi = 2.48 ± 0.06% [28]. 

The 7p cross section for exclusive T(1S) photoproduction was obtained through the rela- 
tion 

where $ is the effective photon flux [18]. The following results were obtained for Q 2 < 1 
GeV 2 : 

a ^T(is) P = 160± 51+48 pb) 60<W<130GeV (6) 
a ^T(is) P = 321 ± 8g +46 4 pbj 13Q < w < 22Q GeV ^ 
a ^r(is) P = 235 ± 47 +so pb) qo<W < 220 GeV (8) 

The number of events, the acceptance, the flux and the cross sections in the different W 
intervals are given in Table [U 

Figure 2 shows the extracted cross section for the two independent W ranges, 60 < W < 
130 GeV and 130 < W < 220 GeV. Also shown are a previous ZEUS result [2] based on 
a partially overlapping data set and the HI result [3]. The measured two cross section 
values were used to calculate 5, resulting in 5 = 1.2 ± 0.8, a value consistent with the 
theoretical expectation [5]. 

The data are compared to several theoretical calculations. Frankfurt, McDermott and 
Strikman (FMS) [5] based their calculation on a two-gluon exchange between the inter- 
acting qq dipoles and the proton, using CTEQ4L parton density functions (PDFs) [29]. 
Ivanov, Krasnikov and Szymanowski (IKS) [7] use a next-to-leading-order (NLO) cal- 
culation in which the prediction for the W dependence of the cross section depends on 
the scale adopted for the hard scattering (presented here for /x =1.3 and 7 GeV). Martin, 
Nockles, Ryskin and Teubner (MNRT) [6] have a NLO calculation using gluon densities 
extracted from HERA data on exclusive J/ip electroproduction. Rybarska, Schafer and 
Szczurek (RSS) [8] use a /^-factorisation approach trying a Gaussian-like and a Coulomb- 
like light-cone wave-function for the vector meson. All the calculations are consistent with 
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the data. In the IKS case, a calculation using an intermediate scale is preferred. In the 
RSS case, the data seem to favour a Gaussian-like wave function. 

9 Summary 

The exclusive photoproduction of T(1S) meson has been studied at HERA with the 
ZEUS detector in the kinematic range 60 < W < 220 GeV, Q 2 < 1 GeV 2 and y/s = 
318 GeV using the muon decay channel. The dependence of the cross section on W 
has been extracted and is in agreement with predictions of several calculations based on 
perturbative QCD. 
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W range (GeV) 


60-130 


130-220 


60-220 


N tot 


159 ±13 


135 ± 12 


294 ± 17 




41 ± 13 


44 ± 12 


85 ± 17 


Nr(is) 


30 ±9 


32 ±9 


62 ± 12 


A 


0.216 


0.230 


0.226 




0.30±0.10i^ 


0.31 ± 0.08±g;?l 


0.60 ±0.12^ 




8.9 ±2.81^ 


9.0 ± 2A±\\ 


17.6 ±3.5^| 




0.055 


0.028 


0.074 


W (GeV) 


100 


180 


140 


^.a^iv.Bi {pb) 


5.5 ± 1.8^ 


10.9 ±3.0^| 


8.1 ± l-6±?;§ 


(j7P ^T(lS)p ( pfc ) 


160 ± 5ltg 


321 ± 88+^ 4 


235 ± 47+^ 



Table 1: The T production cross section for Q 2 < 1 GeV 2 . The first uncertainty 
is statistical, the second systematic. N tot is the total number of events in the signal 
mass region, N S i g is the number of signal events in the signal mass region, iVx(is) 
is the extracted number ofT(lS) signal events, A is the acceptance, Bi is the decay 
branching ratio to muons o/T i; $ is the effective photon flux used to compute the 
7P cross section from the ep cross section at W and a median Ql = 10~ 3 GeV 2 . 
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ZEUS 




Mass |iV (GeV) Mass (GeV) 

Figure 1: Invariant mass distributions of the fi + fi~ pairs in different W regions: 
(a) 60 < W < 220 GeV, (b) 60 < W < 130 GeV and (c) 130 < W < 220 GeV. The 
full dots are the ZEUS data. The hatched and full histograms represent, respectively: 
the GRAPE distribution of the Bethe-Heitler (BH) background and the sum of 
DIFFVM distributions of the T(1S), T(2S) and T(3S) signals. The solid line is 
the sum of the two contributions (see text). 
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Figure 2: The exclusive T(1S) photoproduction cross section as a function ofW. 
The full squares are the ZEUS data from this measurement in the kinematic region 
Q 2 < 1 GeV 2 , and two W ranges, 60 < W < 130 GeV, and 130 < W < 220 GeV. 
The inner bars indicate the statistical uncertainties, the outer bars are the statistical 
and systematic uncertainties added in quadrature. The earlier measurements of 
ZEUS [2] and HI [3] , are also shown. The shaded area denotes predictions of 
NLO MNRT [6]. The long-dashed line is the prediction of the FMS model [5]. 
The dashed-dotted (dotted) line is the prediction of the IKS [7] using a scale of 1.3 
GeV (7 GeV). The small-dashed line (small-dashed three-dots) is the prediction of 
RSS [8], using a Gaussian-like (Coulomb-like) wave function. 
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